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ABSTRACT 

This  paper  focuses  on  evaluating  the  performance  of 
physically  based  distributed  and  conceptual  models,  and 
assesses  their  sensitivity  to  changes  in  the  temporal  and 
spatial  sampling  of  rainfall.  The  Hydrologic  Modeling 
System  (HMS)  was  selected  to  represent  conceptual 
hydrologic  models,  while  MIKE-SHE  and  GSSHA  were 
selected  to  represent  distributed  physically  based  models. 
This  manuscript  presents  results  with  MIKE-SHE,  while 
the  poster  at  the  conference  will  include  results  from 
GSSHA.  The  performance  evaluation  criterion  is  the 
overall  agreement  between  observed  and  predicted 
hydrographs  and  the  models’  ability  to  predict  time  and 
magnitude  of  peak  discharges  and  runoff  volume.  Both 
models  were  carefully  calibrated  and  validated  using 
numerous  storm  events  for  a  21.4  km2  watershed  in 
northern  Mississippi.  The  results  indicated  that  MIKE- 
SHE  captured  the  peak  runoff  discharges  and  total  runoff 
volume  better  than  HMS.  However,  overall,  the 
performance  of  both  models  was  quite  reasonable.  To 
assess  the  models’  requirements  for  rainfall  information, 
an  in-depth  investigation  of  the  impact  of  the  spatial  and 
temporal  sampling  of  rainfall  on  the  prediction  accuracy 
of  each  model  was  conducted.  The  study  showed  that 
MIKE-SHE  was  more  sensitive  to  both  the  spatial  and 
temporal  sampling  of  rainfall  than  HMS. 

1.  INTRODUCTION 

Currently  available  watershed  models  range  from 
simple  conceptual  lumped  models  to  comprehensive 
physically  based  distributed  models.  Conceptual  lumped 
models  use  an  integrated  description  of  parameters 
representing  an  average  value  over  the  entire  catchment. 
A  watershed  can  be  divided  into  a  number  of  sub¬ 
catchments  where  the  hydrologic  parameters  may  vary 
from  one  sub-catchment  to  another.  In  such  case,  lumped 
models  may  be  labeled  as  “semi-distributed.”  They 
remain  non-physically  based,  however,  as  they  use 
synthetic  methods  of  transforming  rainfall  to  runoff. 

Distributed  physically  based  models,  on  the  other 
hand,  can  account  for  spatial  variations  in  input 
parameters  and  state  variables  within  the  catchment.  They 
incorporate  physical  formulations  of  the  different 
hydrologic  processes.  Therefore,  this  class  of  models  has 


the  advantage  of  simulating  complex  hydrologic  systems 
and  utilizing  distributed  field  hydrologic  data.  It  is  also 
recognized  that,  compared  to  conceptual  lumped  models, 
physically  based  distributed  models  are  more  complex  to 
setup,  have  more  stringent  data  requirements,  and  can  be 
subject  to  over-parameterization.  However,  the 
increasing  availability  of  distributed  data  on  rainfall  and 
watershed  properties,  along  with  the  exponential 
improvement  in  computational  resources,  have  increased 
the  interest  of  both  research  and  applied  communities  in 
the  development  and  applications  of  such  models. 

Recent  advances  in  the  development  of  both 
conceptual  and  physically  based  models  have  lead  to  a 
number  of  model  inter-comparison  and  evaluation  studies. 
A  detailed  discussion  of  such  studies  is  given  in  Michaud 
and  Sorooshian  (1994a),  Refsgaard  and  Knudsen  (1996), 
and  Perrin  et  al.,  (2001).  A  review  of  these  comparative 
studies  indicates  that  the  performance  accuracy  of  the  two 
modeling  approaches  varies  widely.  The  comparable 
performance  accuracy  obtained  with  three  models  of 
varying  degrees  of  complexity  lead  Refsgaard  and 
Knudsen  (1996)  to  recommend  the  use  of  conceptual 
models  especially  when  calibration  data  is  available,  and 
to  limit  the  use  of  complex  physically  based  data  for 
ungauged  basins  where  they  are  expected  to  have  a  better 
performance.  As  discussed  by  Refsgaard  and  Knudsen 
(1996),  the  superiority  of  complex  physically  based 
models  over  simpler  conceptual  models  remains  at  the 
hypothesis  level  and  has  not  been  unambiguously 
supported  by  actual  and  sufficient  performance  evaluation 
tests.  Recently,  the  Hydrology  Laboratory  of  the  National 
Weather  Service  (NWS)  office  of  hydrology  has 
conducted  an  extensive  model  inter-comparison  study  to 
assess  the  performance  of  several  physically  based 
models  against  operational  lumped  models  (Smith  et  al., 
2004;  Reed  et  al.,  2004).  The  study  found  that,  in  more 
cases,  the  lumped  model  outperformed  the  distributed 
models.  However,  the  NWS  study  indicated  a  wide  range 
of  accuracies  among  model  results  and  suggested  that 
factors  such  as  model  formulation  and  the  modeler’s  skill 
can  have  bigger  impact  than  the  type  of  the  used  model. 

The  present  study  builds  on  the  continuous  research 
efforts  to  investigate  the  capabilities  and  limitations  of 
conceptual  versus  physically  based  models.  Specifically, 
the  study  evaluates  the  prediction  accuracy  of  two 
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different  rainfall-runoff  models,  conceptual  semi- 
distributed  and  physically  based  distributed  models, 
focusing  on  the  sensitivity  of  both  modeling  approaches 
to  the  quality  and  accuracy  of  the  rainfall  information.  It 
has  been  a  common  belief  by  the  hydrologic  community 
that  rainfall  variability,  both  in  space  and  time,  has  a 
significant  effect  on  the  response  of  hydrologic  systems. 
Therefore,  a  number  of  studies  have  addressed  the  impact 
of  rainfall  sampling  resolution  on  the  prediction  accuracy 
of  hydrologic  models,  and  somewhat  mixed  results  were 
often  reported.  Wilson  et  al.  (1979),  Krajewski  et  al. 
(1991),  Ogden  and  Julien  (1993),  Holman-Dodds  et  al. 
(1999),  and  Michaud  and  Sorooshian  (1994b),  have 
shown  that  the  variability  and  the  sampling  resolution  of 
rainfall  can  have  a  significant  influence  on  the  results  of 
hydrologic  models.  Other  studies  by  Beven  and 
Homberger  (1982)  and  Obled  et  al.  (1994),  however,  have 
shown  that  a  correct  estimation  of  the  total  rainfall 
volume  is  more  important  for  accurate  flow  prediction 
than  providing  the  model  with  detailed  spatial  and 
temporal  patterns. 

The  present  study  has  two  main  objectives:  (1)  To 
compare  the  performance  accuracy  of  two  different 
hydrologic  models,  conceptual  and  physically  based,  and 
(2)  to  investigate  the  relative  sensitivity  of  both  modeling 
approaches  to  the  sampling  resolution  of  rainfall. 


2.  THE  STUDY  SITE 

The  Goodwin  Creek  experimental  watershed  located 
at  the  north  central  part  of  Mississippi  was  selected  for 
this  study.  The  National  Sediment  Laboratory  of  the 
United  States  Department  of  Agriculture  in  Oxford 
Mississippi  has  been  monitoring  the  watershed  since 
1981.  Detailed  information  about  the  watershed  and  the 
available  data  can  be  found  in  Alonso  and  Binger  (2000). 

The  watershed  has  a  fairly  steep  topography  with 
drainage  area  of  21.4  km2.  The  terrain  elevation  in  the 
watershed,  with  reference  to  mean  sea  level,  ranges  from 
71  m  near  the  outlet  to  128  at  the  catchment  divide  with 
an  average  channel  slope  of  0.004.  It  has  a  humid  climate 
(hot  in  summer  and  warm  in  winter),  an  average  annual 
temperature  of  about  65oF,  an  average  annual  rainfall  of 
about  1440.2  mm  (56.7  inches),  and  a  mean  annual  runoff 
of  144.8  mm  (5.7  inches). 

A  network  of  30  gauges  is  used  to  measure 
precipitation  over  the  watershed.  The  watershed  has  been 
divided  into  14  sub-catchments  with  a  flow-recording 
flume  constructed  at  the  outlet  of  each.  Figure  1  shows  a 
map  of  Goodwin  Creek  watershed  along  with  the 
locations  of  the  monitoring  stations.  The  land-use  in  the 
watershed  can  be  described  as  follows:  idle  land  and 
pasture  (60%),  forest  (26%)  and  cultivated  land  (14%), 


while  the  soil  can  be  classified  into  two  main  types:  Silt 
loam  (80%)  and  clay  loam  (20%)  (Blackmarr,  1995). 


Figure  1:  Location  of  the  monitoring  stations  in  the 
Goodwin  Creek  watershed. 


3.  MODEL  CALIBRATION  AND  VALIDATION 

Both  models,  MIKE-SHE  and  HMS,  were  calibrated 
and  validated  using  the  split-sample  test  method  as 
described  by  Klemes  (1986).  Accordingly,  one  set  of  data 
was  used  to  calibrate  the  models  while  another  of  data 
were  reserved  for  model  validation.  During  the 
calibration  procedure,  physical  and  numerical  parameters, 
such  as  loss  and  routing  coefficients  were  adjusted  and 
fine-tuned  to  minimize  the  difference  between  the  model 
results  and  the  field  observations.  However,  parameter 
adjustments  remained  within  the  physically  acceptable 
ranges  based  on  information  available  in  the  literature.  It 
should  be  noted  that  both  models  were  calibrated  using 
the  runoff  discharge  measured  at  the  outlet  of  each  sub¬ 
basin  within  the  watershed. 

Afterwards,  a  second  independent  set  of  data  was 
used  to  validate  the  models.  No  further  adjustments  to  the 
parameters  were  allowed  at  this  stage,  and  the  land  use 
and  other  watershed  characteristics  were  assumed  to 
remain  unchanged. 


4.  EVALUATION  CRITERIA 

The  criteria  used  to  evaluate  the  performance  of  the 
models  are  the  overall  agreement  between  predicted  and 
measured  runoff  discharges,  and  the  models’  ability  to 
predict  time  and  magnitude  of  hydrograph  peaks,  and 
runoff  volume.  The  following  statistical  measures  were 
used  to  quantify  the  performance  accuracy  of  both  models 
during  each  simulation  periods,  and  combined  over  all 
periods: 
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( 1 )  Absolute  Runoff  Volume  Error 
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Where,  Ve  is  the  runoff  volume  error  (m3  or  %),  Vr  is 
the  reference  runoff  volume  (m3),  and  Vp  is  the  predicted 
runoff  volume  (m3). 


(2)  Root  Mean  Square  Error 
Runoff: 
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Where,  <2p  is  the  predicted  Flow  (m3/s),  Qr  is  the 
reference  flow  (m3/s),  Qr  is  the  mean  flow  of  reference 
(m3/s),  i  is  the  hourly  counter,  and  N  is  the  number  of 
discharge  observed-predicted  pairs. 

Rainfall: 
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Where,  Rri  is  the  rainfall  intensity  a  reference 

sampling  frequency,  RAt .  is  the  aggregated  intensity  at 
sampling  frequency  of  At  (30-minute,  1-hr,  2-hr,  etc.). 


(3)  Relative  Peak  Error 


Where,  Pe  is  relative  peak  error  (%),  Pp  is  the 
predicted  Peak  Flow  (m3/s),  Pr  is  the  reference  peak  flow 
(m3/s),  i  is  the  peak  counter,  and  N  is  the  number  of 
Peaks. 

(4)  Error  in  Peak  Time 

N 

Y\t.-t\ 

pi  n\ 

rp  _  1=1  t  r\ 


Where  Te  is  the  error  in  peak  time  (minutes),  Tp 
is  the  predicted  peak  time,  Tr  is  the  reference  peak  time, 
and  N  is  the  number  of  peaks. 

The  reference  quantities  used  in  the  above  statistics 
will  be  clearly  defined  in  the  upcoming  sections. 


5.  CALIBRATION  AND  VALIDATION  RESULTS 

It  is  important  to  state  that  the  raw  precipitation  data 
were  aggregated  to  15  minutes  accumulations  and  used  as 
input  to  both  models.  Figures  2  through  4  and  Table  1 
show  a  summary  of  the  models  performance.  The 
statistics  shown  in  Table  1  are  based  on  the  all  the 
validation  periods  combined  together. 


Figure  2:  Measured  and  simulated  discharge  at  the 
catchment  outlet  (calibration  period). 
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Weighted  Average:  Pe 
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Figure  3:  Measured  and  simulated  discharge  at  the 
catchment  outlet  (validation  periods). 
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Table  1:  Statistical  Summary  of  the  Models  Performance 


Volume  Error 

RMSE 

Relative  Peak  Error 
% 

Peak  Time 
Error 
(min) 

(m') 

% 

m3/s 

% 

Simple  Weighted 

Calibration 

MIKE-SHE 

-28,766 

-1.8 

1.17 

48.3 

3.1 

2.5 

18.0  1 

HMS 

33,213 

2.1 

1.7 

68.4 

5.9 

6.8 

30.2 

Validation 

MIKE-SHE 

726,370 

4.9 

1.9 

63.6 

15.0 

11.9 

30.3 

HMS 

1,116,109 

7.5 

2.1 

69.9 

25.5 

19.1 

24.6 

also  comparable  to  experiments  performed  by  Shah  et  al. 
(1996)  where  a  simple  lumped  and  distributed  models 
performed  well  under  wet  antecedent  moisture  conditions. 


Figure  4:  Measured  and  simulated  peak  runoff  discharges 
during  calibration  and  validation  periods. 

Overall,  MIKE-SHE  reproduced  the  details  of  the 
measured  hydrographs,  while  HMS  was  unable  to  capture 
such  details  due  to  utilizing  synthetic  and  approximate 
runoff  transformation  techniques  coupled  with 
assumptions  of  linearity  and  superposition.  The  statistical 
measures  represented  in  Table  1,  however,  show  that  both 
models  predicted  runoff  volume,  overall  runoff  discharge, 
peak  discharge,  and  timing  of  the  peak  discharge 
reasonably  well.  It  can  also  be  observed  that  the 
distributed  model  performed  better  than  the  conceptual 
model,  especially  in  the  ability  to  predict  peak  runoff 
discharges  (see  Figure  4).  However,  difference  in  the 
performance  between  the  two  models  was  not  drastic.  It 
is  possible  that  the  response  of  both  models  was  similar 
due  to  the  relative  homogeneity  of  the  land  use  and  soil 
properties  of  the  watershed  tested  herein. 

The  performance  exhibited  by  both  models  in  this 
study  is  in  full  agreement  with  other  published  results, 
e.g.  Michaud  and  Sorooshian  (1994a)  and  Refsgaard  and 
Knudsen  (1996).  Both  studies  concluded  that  conceptual 
lumped  and  physically  based  distributed  models  have 
similar  runoff  prediction  accuracy  when  data  is  available 
for  calibration  purposes.  The  results  presented  here  are 


6.  RAINFALL  TEMPORAL  SAMPLING  ANALYSIS 

This  experiment  was  designed  to  study  the  sensitivity 
of  the  models  to  changes  in  the  temporal  sampling  of 
precipitation  data.  The  15-minute  rainfall  data  were 
aggregated  into  30-minutes,  1-hr,  2-hr  and  6-hr 
samplings.  This  aggregation  will  cause  gradual  loss  of 
rainfall  temporal  information  while  conserving  the  total 
volume.  The  predictions  obtained  with  the  15-minute 
precipitation  data  were  considered  to  be  the  reference  to 
which  other  temporal  samplings  are  compared.  The 
simulation  periods  used  to  calibrate  and  validate  the 
models,  were  repeated  for  each  temporal  sampling. 

A  summary  of  the  impact  of  the  temporal  sampling 
on  the  models’  response  is  shown  in  Figures  5  and  Table 
2.  In  order  to  confirm  the  deterioration  pattern  of  the 
predicted  runoff  hydrographs,  the  simulations  were 
repeated  using  a  single  gauge  (gauge#54)  located  at  the 
center  of  the  watershed.  The  results  of  the  test  were 
similar  to  the  all  gauges  temporal  analysis. 


Figure  5:  Error  in  rainfall  and  runoff  due  to  rainfall 
temporal  sampling 
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overall  trend  of  the  response  of  both  models  is  similar, 
MIKE-SHE  was  clearly  and  consistently  more  sensitive  to 
changes  in  the  temporal  sampling  of  rainfall. 


It  is  observed  in  Figure  5  and  Table  2  that  the  runoff 
RMSE,  relative  error  in  peak  magnitude  and  time  vary 
almost  linearly  with  changes  in  the  temporal  rainfall 
sampling.  The  results  also  indicate  that  although  the 


Table  2:  Statistics  Summary  for  the  Temporal  Sampling  Analysis 


|  Rainfall  | 

|  30  min 

2  hr 

6  hr 

WBM 

0.49 

0.60 

e _ 

umi 

199.10 

246.04 

Table  3:  Statistical  Analysis  for  the  Models’  Sensitivity  to  Changes  in  the  Spatial  Sampling  of  the  Rainfall  Data 


Rainfall 


Volume 

mm 

m 

.72 

16.37 

1  -4 

.01 

-6.14 

-20.68 

% 

m 

.06 

1.49 

m 

.37 

-0.56 

-1.88 

Volume  Error 

HS3H 

% 

-0.19 

338,726 

2.25 

i  -96,741 

-0.64 

-201,648 

-1.34 

-505,503 

-3.35 

MIKE-SHE 


HMS 


RMSE 


Error  in  Peak 
Time  (min) 


Relative  Peak  Error 

Error  in  Peak 
Time  (min) 

Simple 

Weighted 

3.24 

2.42 

0.47 

5.26 

4.60 

3.75 

5.68 

5.08 

8.91 

5.66 

5.06 

9.84 

10.10 

9.46 

10.78 
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7.  RAINFALL  SPATIAL  SAMPLING  ANALYSIS 

As  indicated  earlier,  both  models  were  calibrated  and 
validated  using  all  30  rain  gauges.  This  represents  a 
spatial  density  of  about  1 .4  gauges  per  square  kilometer. 
Such  density  was  assumed  to  capture  the  true  spatial 
variation  of  rainfall  over  the  study  area,  as  such,  it  was 
used  as  the  reference  to  which  results  from  other  spatial 
sampling  simulations  will  be  compared  to.  In  all  the 
simulations  presented  in  this  section,  the  15-minute 
rainfall  temporal  sampling  was  used. 

The  number  of  gauges  was  systematically  reduced  to 
20,  10,  5,  2,  and  1  gauge(s)  to  represent  lower  spatial 
sampling  scenarios.  At  every  spatial  sampling  scenario, 
the  selected  rain  gauges  were  as  uniformly  distributed 
over  the  watershed  area  as  possible.  The  calibration  and 
validation  simulation  periods  were  repeated  for  the 
various  spatial  sampling  experiments.  It  should  be 
emphasized  that  the  error  introduced  into  the  rainfall  data 
affects  both  its  spatial  variability  and  total  volume. 


8.  COMBINED  SPATIAL-TEMPORAL  SAMPLING 
ANALYSIS 

In  this  set  of  experiments,  both  temporal  and  spatial 
samplings  were  allowed  to  vary,  while  experiment  with 
the  highest  temporal  sampling  of  15  minutes  and  spatial 
sampling  of  30  gauges  was  used  as  the  overall  reference. 
For  each  temporal  sampling  of  1 5-minutes,  30-minutes,  1- 
hour,  2-hours,  and  6-hours,  a  simulation  with  1,  2,  5,  10, 
20  and  30  gauge(s)  was  performed.  The  main  objective 
of  this  set  of  experiments  is  to  observe  the  response  of 
both  models  to  combined  changes  in  the  temporal  and 
spatial  sampling  of  rainfall.  Table  4,  and  Figures  7  and  8 
show  the  RMSE  of  the  runoff  prediction  as  function  of 
different  spatial  and  temporal  sampling.  The  relative 
error  of  peak  magnitude  and  time  showed  similar  patterns 
to  those  of  the  runoff  RMSE. 


Table  4:  RMSE  (%)  for  Combined  Spatial-Temporal 
Analysis. 


The  statistical  analysis  for  these  experiments  is 
summarized  in  Table  4  and  Figure  6.  The  results  show 
that,  generally  speaking,  the  models’  performance 
deteriorates  as  the  density  of  rain  gauges  in  the  watershed 
is  decreased.  Moreover,  the  figures  also  show  that  the 
response  of  both  models  to  changes  in  the  rain  gauges 
density  was  somewhat  similar  with  MIKE-SHE  showing 
higher  sensitivity  especially  at  the  lower  end  of  the  spatial 
sampling.  For  example,  increasing  the  number  of  rain 
gauges  from  2  to  5  has  lead  to  a  reduction  in  the  RMSE  of 
MIKE-SHE  from  32%  to  22%,  while  the  RMSE  of  HMS 
decreased  from  22.4%  to  19.6%  only. 
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Figure  6:  Error  in  rainfall  and  runoff  due  to  changes  in 
spatial  sampling  of  rainfall. 
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Figure  7:  Runoff  error  distribution  due  to  rainfall 
temporal-spatial  resolution  (MIKE-SHE). 
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Rainfall  Resolution  (min) 

Figure  8:  Runoff  error  distribution  due  to  rainfall 
temporal-spatial  resolution  (HMS). 

Figures  7  and  8  show  that  the  response  of  the  models 
to  changes  in  the  rain  gauges  density  is  a  function  of  the 
temporal  sampling  at  which  these  gauges  are  set.  For 
example,  the  impact  of  decreasing  the  density  of  rain 
gauges  is  minor  if  they  are  set  to  temporally  sample  at  6- 
hour  intervals.  On  the  other  hand,  and  as  clearly  shown  in 
Figures  7  and  8,  if  the  rain  gauges  are  set  to  temporally 
sample  at  short  intervals  (1-hour  or  less),  there  is  a 
significant  deterioration  in  the  models’  performance  as 
the  density  of  the  rain  gauges  decrease.  It  can  also  be 
observed  that  5  gauges  set  to  temporally  sample  at  15- 
minutes  deliver  an  equivalent  model  performance  to  a  10- 
gauge  network  set  to  temporally  sample  at  1-hour 
intervals.  In  other  words,  one  might  compensate  for  the 
loss  of  rainfall  spatial  information  by  increasing  the 
temporal  sampling. 


9.  CONCLUSIONS  AND  CLOSING  REMARKS 

Physically  based  distributed  model  (represented  by 
MIKE-SHE)  and  conceptual  semi-distributed  model 
(represented  by  HEC-HMS)  performances  were  compared 
in  this  study.  Both  models  were  setup,  calibrated,  and 
validated  for  the  Goodwin  Creek  watershed  in  northern 
Mississippi.  The  impact  of  the  temporal  and  spatial 
sampling  of  rainfall  on  the  performance  of  both  models 
was  also  investigated. 

The  simulations  performed  in  this  study  focus  on  two 
components  of  errors.  The  first  component  is  the  error 
related  to  factors  such  as  model  structure,  formulation, 
lack  of  sufficient  data,  etc.  This  error  characterizes  the 
overall  model  performance  and  was  quantified  herein  by 
comparing  the  observed  measurements  to  model 
simulations  obtained  with  the  15  minutes  temporal 
sampling  and  all  30  rainfall  gauges.  The  second 
component  of  error  characterizes  the  sensitivity  of  model 
performance  to  the  temporal  and  spatial  rainfall  sampling. 


In  evaluating  their  performance,  both  models 
predicted  runoff  volume,  overall  runoff  discharge,  peak 
magnitude  and  time  with  a  reasonable  accuracy.  The 
error  in  predicting  runoff  volume  was  less  than  8%, 
RMSE  was  less  than  70%,  and  the  error  of  peak 
magnitude  was  less  than  26%.  The  statistical  measures 
showed  that  MIKE-SHE  predicted  the  runoff  volume, 
overall  runoff  discharge,  and  peak  magnitude  better  than 
HMS.  On  the  other  hand,  HMS  outperformed  MIKE- 
SHE  in  predicting  the  peak  timing.  However,  the  results 
of  both  models  were  overall  quite  comparable.  This 
realization  is  in  agreement  with  the  works  of  Michaud  and 
Sorooshian  (1994)  and  Refsgaard  and  Knudsen  (1996). 
Both  studies  concluded  that  conceptual  and  physically 
based  models  have  similar  runoff  prediction  accuracy 
when  data  is  available  for  calibration.  The  impact  of 
employing  linear  conceptual  methodologies  in  HMS  was 
evident  in  its  inability  to  capture  double-peak  rainstorms, 
and  its  sensitivity  to  any  changes  to  the  rainfall  volume. 
It  is  noteworthy  that  identifying  and  using  seasonally 
variable  calibration  parameters  was  not  the  focus  of  this 
study.  Therefore,  all  simulation  periods  used  in  the 
calibration  and  validation  were  limited  to  the  wet  non¬ 
growing  season  of  January  through  May.  Accordingly,  it 
was  possible  to  use  a  single  set  of  calibration  parameters. 

In  reference  to  the  sensitivity  of  the  models  to  the 
spatial  and  temporal  rainfall  sampling,  the  study  showed 
that  errors  introduced  by  coarse  sampling  scenarios  can  be 
significant.  For  example,  the  second  component  of 
prediction  errors  was  in  the  order  of  V2  to  of  the  first 
component  of  the  error  as  a  result  of  reducing  the 
temporal  sampling  from  1 5  minutes  to  2  hours.  Similarly, 
reducing  the  number  of  gauges  from  30  to  2  resulted  in  a 
second  component  of  prediction  errors  in  the  order  of  to 
V2  of  the  first  component  of  the  error.  Overall,  for  this 
particular  watershed  size,  increasing  the  rain  gauge 
density  from  1  to  2  resulted  in  the  most  significant 
improvement  for  both  models.  Similarly,  a  temporal 
sampling  frequency  beyond  1  hour  showed  significant 
deterioration  in  the  quality  of  the  runoff  prediction. 

This  study  also  showed  that  MIKE-SHE  was  more 
sensitive  to  the  rainfall  temporal  and  spatial  sampling  than 
HMS.  Such  sensitivity  was  more  pronounced  and 
persistent  especially  when  the  spatial  sampling  was 
significantly  lowered.  The  sensitivity  of  MIKE-SHE  can 
be  attributed  to  its  inherent  dependency  on  the  spatial 
distribution  of  input  data,  and  the  physically  based 
methodologies  employed  to  model  the  various 
components  of  the  hydrologic  cycle.  This  observation 
emphasizes  the  need  for  detailed  rainfall  information  to 
obtain  accurate  runoff  prediction  using  distributed 
physically  based  models. 
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The  combined  spatial-temporal  sampling  experiment 
showed  that  increasing  the  temporal  sampling 
compensates,  at  least  partially,  for  the  loss  of  rainfall 
spatial  information.  It  also  showed  that  under  poor  spatial 
sampling  conditions,  the  gain  in  model  performance  by 
increasing  the  temporal  sampling  frequency  becomes 
negligible. 

It  should  be  emphasized  that  the  conclusions  and 
findings  of  this  study  are  conditioned  on  the  hydrologic 
and  meteorologic  characteristics  of  the  selected 
watershed.  Expanding  the  experiments  performed  herein 
to  include  watersheds  of  various  sizes,  and  rainfall- 
runnoff  response  characteristics  would  further  evaluate 
the  performance  of  each  model  and  its  sensitivity  to 
temporal  and  spatial  sampling.  Efforts  are  underway  to 
evaluate  the  performance  of  conceptual  and  physically 
based  models  on  coastal  low-gradient  watersheds  with 
high  variable  tropical  rainfall  regimes. 
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